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Abstract 


The  steady  matian  af  a  planing  surface  of  moderate  aspect  ratio  at  small  angles 
af  attack  is  considered.  Linearized  theory  is  used  with  a  square-root  tyj>e  of  pressure 
singularity  representing  the  flow  near  the  leading  edge.  An  asymptotic  solution  for 
the  pressure  distribution  on  the  planing  surface  at  large  Froude  number  (or  small  ^  , 
the  inverse  af  the  Froude  number)  is  sought.  The  lowest  order  term  of  the  pressure 
distribution,  obtained  by  setting  ^  equal  to  zera,  is  found  to  be  the  same  as  the 
pressure  distribution  on  the  lower  side  of  the  corresponding  thin  wing.  Higher  order 

terms  in  ^  are  obtained  by  an  iteration  process.  Explicit  solutions  are  obtained  ta 

2 


order  ^3  for  rectangular  planforms.  Numerical  results  are  calculated  for  rectangular 
flat  piate  planing  surfaces  of  aspect  ratios  from  0.5  to  2.0.  It  is  foi  .d  that  for  large 
aspect  ratios  the  lift  coefficient  is  reduced  by  the  gravity  effect  and  for  small  aspect 
ratios  it  is  increased,  fhe  dividing  aspect  ratio  being  about  1.5.  The  results  compare 
reasonably  well  with  experimental  data. 
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I.  fotroduction 

When  a  surface  craft  moves  at  iow  speed  through  water,  the  lift,  which  supports 
the  craft  on  the  water  surface,  is  supplied  mainly  by  buoyuncy.  If  the  speed  of  the 
craft  is  increased  so  that  the  water  surface  separates  smoothly  from  the  trailing  edge  of 
the  craft,  the  craft  is  said  to  be  planing,  or  gliding,  on  the  water  surface.  During 
planing  motion,  the  lift  is  mainly  supplied  by  hydrodynamic  forces. 

An  important  feature  of  a  planing  motion  is  the  splash  phenomenon,  which  is 
a  spray  sheet  thrown  out  ahead  and  sideways  of  the  planing  craft.  If  the  angle  of 
attack,  which  may  he  defined  as  a  characteristic  angle  between  the  wetted  surface  of 
the  planing  craft  and  the  undisturbed  water  surface,  is  small,  the  thickness  of  the 
splash  is  expected  to  be  thin.  Green  (1935,1936)  made  non  ~  linear  studies  of  two- 
dimensional  flat  plate  planing  at  an  angle  of  attack  of  with  an  infinite  Froude  number. 
The  Froude  number  Fr  is  defined  as  the  ratio  of  the  inertia  effect  to  the  gravity  effect, 
or  Fr  =  ,  where  U  is  the  speed  of  the  planing  craft,  g  the  gravitational  accelera¬ 

tion  and  1  c  characteristic  length.  In  his  results,  the  thickness  of  the  splash  was  found 
to  be  proportional  to  d  for  d  small.  Wagner  (1932)  studied  both  vwo-dimensional 
and  three-dimensional  planing  problems  at  infinite  Fr.  In  his  linearized  formulation, 
the  governing  equations  were  shown  to  be  the  same  as  those  found  in  flows  past  thin 
wings,  except  that  in  planing  problems  the  fluid  under  consideration  is  in  contact 
with  the  wing  on  the  lower  side  only.  It  is  well  known  that  in  thin  airfoil  theory  the 
pressure  has  a  square-root  singularity  at  the  leading  edge  of  the  foil.  Based  on  a  local 
flow  study,  Wagner  showed  that  to  represent  the  splash  in  planing  problems,  the  same 
type  of  pressure  singularity  should  be  used. 
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'Alien  the  effect  of  gravity  is  considered,  two-dimensional  planing  surfaces  of 
various  shapes  have  been  studied  by  many  authors  (cf.  Wehausen  &  Laitone  (I960)). 
Recently,  Mcruo  (1967)  considered  three-dimensional  planing  surfaces  of  large  and  small 
aspect  ratios;  however,  for  small  aspect  ratio,  his  method  requires  Fr  to  be  very  large 
and  is  not  applicable  to  a  rectangular  planing  surface.  In  the  above  solutions  it  was 
assumed  that,  as  in  the  case  of  infinite  Fr,  the  splash  was  a  second  order  quantity  in 
angle  of  attack  and  might  therefore  be  neglected  in  the  formulation  of  the  linearized 
theory  and  that  the  pressure  has  a  square-root  singularity  at  the  leading  edge. 

In  this  paper,  we  consider  a  steady,  three-dimensional  potential  flow  pasta 
planing  surface  of  moderate  aspect  ratio  at  a  large  Froude  number.  We  assume  that 
the  angle  of  attack  is  smal*  so  that  the  problem  may  be  linearized.  In  formulating  the 
problem,  we  represent  the  planing  surface  by  an  unknown  pressure  distribution  over 
the  part  of  the  water  surface  directly  underneath  it.  The  geometric  configuration  of  the 
splash  will  be  neglected  and  the  pressure  is  assumed  to  have  a  square-root  type  of 
singularity  at  the  leading  edge  of  the  planing  surface.  The  perturbation  potential  due 
to  this  pressure  distribution  is  expressed  in  the  form  given  by  Peters  (1949).  It  involves 
a  quadruple  integral  with  the  integrand  linearly  proportional  to  the  unknown  pressure. 
The  kinematic  boundary  condition  on  the  planing  surface  will  lead  to  a  linear  integral 
equation  for  the  unknown  pressure  distribution.  To  facilitate  the  solution  of  the  integral 
equation,  we  expand  its  kernel  for  points  nearby  the  planing  surface,  asymptotically 
for  Fr  — co  ,  In  this  expansion  the  unknown  pressure  distribution  is  regarded 
as  if  it  were  independent  of  Fr,  Then  in  a  similar  fashion  an  iteration  process  will 
yield  successive  terms  in  the  pressure  expansion.  When  Fr  is  set  to  infinity,  only  the 
lowest  order  term  in  the  expansion  remains.  This  term  corresponds  to  the  downwash 


3 


integral  equation  in  wing  theory,  except  that  the  pressure  is  equal  to  one-half  of  the 
loading  on  the  corresr^r.c-  iny.  The  downwosh  integral  equation  is  solved  using  a 
method  similar  to  that  of  Watkins,  Woolstonand  Cunningham  (1959).  From  the  solution 
of  the  downwash  integral  equation,  the  lowest  order  term  of  the  pressure  expansion  in  Fr 
is  obtained  in  terms  of  the  given  "downwash".  At  each  stage  of  the  iteration  process  the 
same  dowr.wash  integral  equation  has  to  be  solved.  The  solution  of  the  downwash  inte¬ 
gral  equation  provides  each  new  term  in  the  expansion  of  the  pressure  in  terms  of  the  pre¬ 
viously  obtained  terms.  Our  solution  is  carried  out  up  to  terms  of  order  Fr“2.  The 
above  iteration  scheme  has  been  advanced  by  Cumberbatch  (1958)  in  solving  twa-dimen- 
sional  planing  problems  at  high  Froude  number. 

In  this  paper  numerical  results  are  given  for  rectangular  flat  plate  planing  sur¬ 
faces  having  aspect  ratios  from  0„5  to  2.0. 

2.  Derivation  of  the  integral  Equation 

Consider  a  pianing  surface  of  moderate  aspect  ratio  travelling  at  a  constant 
velocity  U  over  a  water  surface  of  infinite  extent.  The  angle  of  attack  is  assumed  to 
be  small,  so  that  linearized  theory  may  be  adopted.  The  water  is  considered  to  be  in- 
viscid,  incompressible  and  of  infinite  depth. 

We  choose  a  set  of  Cartesian  coordinate  axes  x-y-z  fixed  to  the  planing  surface. 
The  x-y  plune  is  assumed  to  coincide  with  the  undisturbed  water  surface,  the  x-axis 
pointing  in  the  direction  opposite  to  the  velocity  U  and  the  z-axis  in  the  direction  op¬ 
posite  to  the  gravitational  acceleration  g.  In  this  frame  of  reference,  the  fluid  at  in¬ 
finity  appears  to  have  a  uniform  velocity  U  in  the  x -direction  and  the  motion  becomes 
steady.  We  choose  U  as  our  characteristic  velocity.  For  planing  problems  the  proper 
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characteristic  leng  h  should  be  measured  in  the  flow  direction.  For  a  three-dimensional 
problem  anothe-  characteristic  length  is  the  span  width.  Therefore  two  Froude  numbers 
can  be  def.ned.  For  convenience  in  numerico!  evaluation  we  choose  the  semi-span 
width  as  the  characteristic  length.  Since  we  are  dealing  with  moderate  aspect  ratios 
only,  the  choice  is  not  important.  Based  on  these  characteristics  quantities,  the  velocity 
potential  (x,y  .z)  can  be  written  as 

lU.'J.D  =  X  +  ,  (1) 

where  x  denotes  the  velocity  potentiol  corresponding  to  the  uniform  velocity  when  viewed 
from  the  x-y-z  system  and  (x,y,z)  is  the  perturbation  potential  due  to  the  presence 
of  the  plan'ng  surface.  It  is  obvious  that  4>  satisfies  the  Laplace  equation 

A4>  =  O  .  (2) 

Let  us  denote  the  area  projected  by  the  planing  surface  cn  the  x-y  plone  by  A.  If  we 
represent  the  planing  surface  by  an  unknown  pressure  distribution  on  A,  the  linearized 
boundory  condition  of  may  be  expressed  as  (see  vVehousen  &  Lairone  (I960)), 

lx  *  *  +  2  p  |f-  ”  0  i  on  1*0,  (3) 

where  the  non-dimensional  small  parameter  j3  is  defined  by 

/3  =  ?i/u"  =  l/Pr  <4> 

and  where  jp  =jp(x.,  is  the  non-dimensional  pressure.  We  take  the  ambient  pres¬ 
sure  to  be  zero  ond  define  ^  as  the  ratio  of  the  pressure  to  the  dynomic  head  j>V/ 
where  P  is  the  density  of  the  water. 
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From  the  definition  of  p  it  is  clear  that  on  the  free  surface,  or  on  the  part  of  the  x-y 
plane  outside  A, 


P=°> 


(5) 


and  on  A,  as  stated  before,  p  is  unknown.  The  other  boundary  conditions  on  4> 
may  be  stated  as 


V4>  =  O 


r  X  m  -  t>0 

(6) 

at  )  ^  =  ±  OG 

(7) 

'  i  *  -  oo 

(8) 

Condition  (6)  assures  us  that  no  gravity  wave  will  propagate  upstream. 

The  potential  which  satisfies  the  boundary  conditions  (3),  (5),  (6),  (7)  and  (8) 
is  (see  Peters  (1949)) 


=  ?.K j)  ^  c vsfiypjJ' 

A  o 

X[&Z  costl  -C _  r  f - ,  n  , 


-N(x-|3 


>  (9) 


where  sgn  (x-  §  )  is  the  sign  function  of  x-  £  ,  H(x-  £  )  's  the  Heaviside  step  function 

and 


(10) 
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If  Hie  prof'le  of  the  planing  surface  is  expressed  by 


Z  =  iu.'j)  , 


and  if  we  denote 


W(x,4)  -  li. 
} 


the  linearized  kinematic  boundary  cond'tion  on  the  planing  surface  requires  that 


W(*,M  )  =  I  'm  -|— 
J  l-*0  «  *- 


^ar  { x,  y  j  on  A  .  (13) 


Equations  (9)  and  (13)  will  give  us  an  integral  equation  for  the  unknown  pressure  distri¬ 
bution  on  A.  When  we  differentiate  (9)  with  respect  to  z  a  term  containing  sin  tz  in 
the  t  integration  is  produced.  It  can  be  shown  that  this  term  gives  no  contribution  as  z 
tends  to  zero.  Therefore,  the  integral  equation  for  the  unknown  pressure  distribution  on 


A  becomes 


where 


W<‘'V  =  U  i  . 

=-4z'sTl'l>\ 

0  ^0 

r^^ivr  cc*  >  <^r 


(15) 
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!n  obtaining  (15),  (10)  has  been  used. 


3.  Large  Froude  Number  Expansion 

In  this  section  we  shall  exaand  asymptotically  the  kernel  function  K(x-  |  ,  y-  n 
z;  ^  )  for  large  Fr,  or  small  ,  when  (x,y)  is  on  A  and  z  is  small .  The  method  of 
ex.pcnsion  depends  on  the  aspect  ratio  being  of  order  1  or  smaller. 

Let  us  denote  the  double  integral  occuring  in  K,  shown  in  (15)  by 

r  F  i  rf”  ,, , 

I, ~  e«PH  z''  |xi1Jg-  (,6) 


A  change  of  variables 


X  =  -Rc£>s©  ,  -J  = 


A  0 


transforms  (16)  into 


.'till  rc0 


r  r  0  3 

I,  -\  cos2edei  ~oTz — r  ^(■^li^^e)cos(42cos&)e<p(-iliijcl4 

J  J  ii  +^cose  J 


2  t* 

i 

3 


-  yRe\  cosed© 


where 


2  » 

Z\  )  •/*“* Wpit4} 


>  (17) 


-  I X I  +  i  (^tslne  +  2  cos  © )  , 


*4^  ~ I  *  |  +  i  ( s> t"\(5  —  loose) 


(18) 


The  integrations  with  respect  to  k  in  (17)  can  aii  be  carried  oat  and  we  have 
■■Till  2 

I*  ~  CCS  9 J©  r  |  “~2  -^cose  cx>s(^4,.£g,se)ci(^'^cei&) 

0  *  *  ^ 

-  cc*e,)5i  (jgl^  C&se)j^  ,  ^ 

where  ci(x)  and  si(x)  are  the  cosine  and  sine  integrals  respectively  and  ore  defined  as 


Cc/Xj  =  =  -  T-lms  +  £  (-0 


p*41  x2n 


2n(2n)» 


s c(»)  x  - 


CC  2«  *>l 

Tt +  p  r_jT_ 

z  '  L  1  U  (2r.-t)an-t;| 
n=t 


In  (20),  T  <s  Euler's  constant.  Under  the  assumption  that  the  aspect  ratio  is  not  large 
x[»  ,  "4^  =  0(1)  and  therefore  Jj  may  be  expanded  asymptotically  for  ^  small. 

The  expansion  valid  up  to  the  order  of  fi*  is 

( 

rTcfl  2  j 

It  =  \  oosZ8cleV  f  -qjj  -jfcoso^-T  -  ini^^cose)  ^C(|S4j 


a*!  o 


l[*'Z  i.  ,  f  X*-(H5’»neHa?se)Z  xx-(4<ine  -la>se) 

2  J  cosed  [  (xi  +  (t|i'n0  +  "2oo58)2]Z  !**+( jsine  -  icotTej1]1 

0 

+  B  cos  ©[iT  +  2  Itncosfi  +  Z t  j  In  (XZ  +  [cjs’me 


+  la^aJJix't-tj-sine  -2CO5  0]1)]  +  C  (p3)  J  •  (22) 

In  the  last  expression  of  (22),  the  first  two  terms  inside  the  curly  brackets  are  independent 
fi  ,  and  hence,  gravity  free;  the  remaining  terms  are  dependent  on  ^ 


.  Since 
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we  expect  to  obtain  the  downwash  integral  equation  in  wing  theory,  er.cepf  for  a  factor 
of  1/2 -  from  the  integral  equation  (14),  and  since  when  x  c  £  the  kernel  of  the  integra 
equation  { 14)  is  proportional  to  Jj ,  we  expect  to  recover  the  kernel  of  the  downwash 
integral  equation  for  x  5  from  5'ovity  free  terms  in  (22).  This  means  the*  if  we 
set  z  =  0  in  (22)  we  should  expect  a  singularity  in  y  of  the  form  1/y*  to  appear  after  the 
&  integration.  The  meaning  of  this  type  of  singularity  is  well  known  (Thwaites  (I960)). 

Since  the  remaining  ^  -<■  dependent  terms  in  (22)  will  not  produce  a  singularity  stronger 
than  that  of  the  gravity  free  terms  we  may,  therefore,  set  z  =  0  first  and  then  integrate 
with  respect  to  6.  This  yields 


h1'  + 

T ( 1*7  +  /I  +  0(^}  • 


(23) 


Now,  let  us  denote  the  single  integral  term  occurring  in  (15)  by 

^ (0  ? 


T  -  r 

z 


J 


(24) 


\  c  ^  cos  %)eY?(//'Z)^y  > 
where  is  given  by  (10).  We  note  that  in  (24)  we  have  to  consider  only  x  >  0  as  can 
be  seen  clearly  from  (15).  Since  Jq,  is  an  even  function  of  y  we  need  consider  only  y  >  0. 

The  integral  may  be  regarded  as  representing  part  of  the  downwash  due  to  a 
concentrated  pressure  point  moving  on  the  surface.  It  'nerefore  contains  a  singularity 
of  high  order  aiong  the  track  line  y  =  ^  for  2  =  0,  However,  as  pointed  out  by  Lamb 
(1934),  the  singularity  is  due  'o  the  artificial  nature  of  a  concentrated  pressum  point 
and  disappears  for  a  diffused  pressure.  Therefore,  when  necessary,  the  -  integra- 
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tion  will  niwoys  be  carried  out  first  to  remove  the  singularity  in  (y  -  ^  )  thot  would 
otherwise  appear  as  z  goes  to  zero.  In  the  following  expansion  procedure  for^  we  shall 
assume  z  to  be  different  from  zero,  when  necessary,  in  order  to  avoid  the  oppeoronce 
of  the  singularity  in  y. 

By  using 


(25) 


equation  (24)  may  be  written  os 

-tfi 

2  i 

I,  = 


Let  us  odd 


i  ^  chyU  .  (26) 


to  (26't  to  obtoi.i 


=  j  [^S2etyt*  +  i|8  (ydytziyu  -  x  ci^uj]^. 


(28) 


The  path  of  integration  in  (28)  may  be  changed  to  C  which  runs  from  -  60  +  c  to 


<?T  , 

c<i  +  — •  ^  porallel  to  the  Re  p-axis.  Along  C 


(29) 


where  t  is  purely  real.  The  expansion  of  exp(-i  ^  xchu)  into  power  series  of  (-i ^  xchu) 
and  the  substitution  of  (29)  into  (28)  give 


E 

f 


z  C°°  >'M 

=  /6^p(4p)  Z(^AsU  -t'cUf* 


C<f  f”  "§  (ycbli  -<-Zsl\2 i)]Jt  . 


(30/ 


By  the  binomial  expansion  of  (sht  -  icht)n+^  and  by  separation  into  real  and  imaginary 
Darts  we  can  write 


oo 


K+4 


\  r  (-o‘ 


rti-4  |  n  ->-4-w  _  r“i 

SK  t  ch  i  ce>n(jS2Shzi)  4-  ^  (-  I J  Z  6 


-oo  -m-oy2^--- 
£tH4  w-l 


rt+4 


skn'4'r,rt3'.,(i- '  ■ 


n  ( Ip l skz-t ) j  e *p  (-  J jS^l ^ *.4  )  A  . 


(31) 


It  can  be  seen  that  when  n  is  odd  the  integrand  in  (31)  is  an  odd  function  of  t.  Therefore 
the  only  contribution  is  from  even  values  of  n,  in  which  case  we  can  write 

n+4  jvi-4-w 


I  n-*,4-nt  .  m  I  2  - 7 - 

sk  ick  t  =({)  c  ck  2.1  -  I  )  (ck U 


21 

2. 


rv+4 

I  \  2- 


n+5-w 


i )  y r\  zven 

(32) 


I  2.  2. 

=  (t)  (cki-IJ  (ckz4  +  l)  skz-t  no  otU  . 


If  we  set  n  =  2  V  with  m  =  2p  for  m  even  and  m  =  2p  +  1  for  m  odd  we  get 
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I*  -f  (TTTl^  *  {  £  (->/“ C^dut-O^ 

I  -  3  O  ' 


/(  T  c  ii  / 

co*(-£fi2sk l[)  4-  £  (-J j>  C^;t  (ctai-U  (cKl-4^0  • 

/ 

sVu4  i‘^(Y|6Jsk2t)J  e*p  (-  {jSjeiu-fcjd-l  .  (33) 


In  (33),  we  write 


(cut  t|), 


(34) 


and  noting  that  the  integral  is  uniformly  convergent  for  y  ^  0  we  get 


cos 


(35) 


where  accounts  for  the  factor  sh2t.  The  integral  occuring  in  (35)  can  be 

transformed  into 


I 


2  Re 


(36) 


1 
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where  §  =  tan  \z/y).  For  z  ^  Owe  have  — rr/2  ^  S  0  since  we  need  con¬ 
sider  only  y  >  0.  By  the  change  of  variable 


V  =  it  +  i-S 


(37) 


X  becomes 


- o  iS 


j  «p(-{/|JP7FcU;Jv.  (35) 

»'$  o  ' 

Since  the  first  integral  in  (38)  is  purely  imaginary  and  since  —  Tc/2  ^  O  / 

we  have 


J  =  ^  e* p  (-  jj$ z+  2*  cUm)^v . 


(39) 


The  above  integial  can  be  integrated  (Watson  (1944))  to  give 


(40) 


where  Kc  is  the  Bessel  function  of  the  second  kind  with  imaginary  argument,  which,  for 


j, 3  small  and  J  y^  +  =  0(1),  can  be  expanded 


as 


I  >  -TV>=  O 


2>i 


-I- 


y  i) , 

'  ■  (w! ) 


■WzC 


(41) 


L 


4>{™+d  =  ~ 


4U)  =  -  v* 


ond  iT  is  Euler's  constant.  With  the  aid  of  (40),  I2  given  in  (35)  can  now  be  written 
as 

xi  (|  +  4;  ^ 

- £f+  Z  (-'“clil* i^*f) 

/*«  / 

Ko  (zfi]f+?)  . 

I  (43) 

Using  (41)  and  (42),  we  may  expand  (43)  for  jg  small  as 

T  -  L  c>i\\  ZUll - ofi.ly^—) ^ — 

-1  -  v  ip-'j  P  '*  +  2  *  vy1'  'j'+i2- 

-  If  tln4' +r  +  1 _  + 


+  (x  Ip  +  7e  x4^s)ln  +0C&1) 


18  -^4 


The  substitution  of  (44)  and  (23)  into  (14)  gives 
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Wl,,,3>  *  t1  (f +  h  L;r»S  jH(’'-5>t 

A  A 

1+ Y^-5)  ay-]  Oj-^Ai1  PL*>Vh\  +  +r+-£ 

‘ln4)^  P<MM  -SW'-*) 

+  ^H(X-S)[|h(3-?)2-  -|.]p«.t><s^  +Uo^HOc-ij[z„-^ 


+  (X-5 J|4~  ] I n [ fy ^ i* ] P(|y ^ J-7 J  40(|S5)  ,  (45) 


where 


QU.  J)  =  U  (Ixl  +  Wx)  I  _ 


(46) 


Now,  we  interchange  the  order  of  integration  and  differentiation  with  respect  to 
y  in  (45)  and  then  let  z  tend  to  zero.  In  doing  so  we  must  evaluate 


U  '  f‘U>  7 Z^TTT^l, 


(47) 


J„j  •*.  (j-y 
where  the  integration  is  carried  out  over  the  entire  span,  and  it  is  understood  that  z 
tends  to  zero  through  negative  values.  The  result  is 


-TCpfS^).  (48) 

This  procedure  is  carried  out  to  allow  the  distributed  nature  of  the  pressure  ro  overcome 
the  high  order  singularities  in  y  that  were  mentioned  previously  in  connection  withl2* 
This  gives 


*6 


where  G(  x-  £  ,  y  -  £  )  is  given  in  (46),  the  lower  limit  L.E.  indicates  thot  the 

integration  storts  from  the  leading  edge  of  the  planing  surfoce,  ond  Ax  is  the  part  of  A 
bounded  between  the  leoding  edge  and  the  l?ne  ^  =  x. 


4.  Application  to  a  Planing  Surface  of  Rectangular  Plonform 

In  this  section,  we  sholl  consider  a  rectongular  planing  surface  having  the 
ieading  edge  at  x  =  -b  and  the  trailing  edge  at  x  =  b.  The  integral  equation  shown  in 
(49)  will  be  solved  by  an  interation  process.  The  unknown  pressure  distribution  on  A  is 
expended  into 
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where  i\)  is  the  Birnbaum  expansion  (1923)  derived  from  thin  airfoil  theory. 
It  is  defined  as 


=  sinnfl 


when  n  =  0 


when  n  5- 1 


(51) 


where 

G  =cos"l(-§/b3  • 


(52) 


However,  in  actual  calculation,  the  series  given  in  (50)  will  be  truncated. 

The  iteration  process  starts  by  assuming  ^3=0.  Equation  (49)  reduces  to 


- 


1 


(53) 


Except  for  a  factor  of  1/2,  this  equation  is  the  downwash  integral  equation  for  a  thin 
wing.  Expanding  p(  ^  ^  )  in  the  form  of  (50),  (53)  can  be  solved  numerically.  The 

method  used  is  similar  to  that  of  Watkins,  Woolston  and  Cunningham  (1959).  The 
solution  to  (53)  gives  the  first  term,  p^(  ^  ^  ),  in  the  expansion  of  the  pressure 

as 


M  N 

w.v-Z  L 


_  icr)  y>\ _ _  • 


=  P  W  -  4> 


(54) 


where  M  and  N  are  constants.  Using  the  solution  for  pQ  (  ^  ,  ),  equation  (49), 

when  approximated  to  order  A  may  be  written  as 


!8 


Tp[^P.<5.H >J5  +  (,t'5)lP*(5'3)£!5] 

=  '■'*  jT* -i?ip 


To  solve  (55),  we  write 


f<i.p  =f.(j,z;^pl(|,7) 


The  substitution  of  (56)  into  (55)  and  the  use  of  (53)  and  (54)  give 


-b  J  4 


xi=]  WL12J  , 

•t+n-i)'1  ,yrn  L  ’ 


Jtx-sy’+oj-^/11  <jj-  C'  i3/' 

A 

where  the  integrations  on  the  left-hand  side  con  all  be  carried  out  analytically.  This 

integral  equation  for  p,  ($,  r£  )  is  of  the  same  form  as  (53),  However,  the  second 

2  ~3  /2 

term  on  the  left-hand  side  of  (57)  contains  singularities  of  the  form  (1  -  y  )  '  as 
y  — *.  £  1,  at  the  tips  of  the  planing  surface.  This  singular  behaviour  of  the  in¬ 

duced  downwash  mo;'  be  due  to  the  linearization  of  the  problem.  We  assume  that  away 
from  the  tips  the  error  due  to  linearization  is  smaller.  Also,  the  numerical  methods 
used  to  solve  downwash  integral  equations  of  the  form  (53)  do  not  use  the  values  of 
the  downwash  W (x,y)  at  the  tips.  Furthermore  since  the  pressure  vanishes  at  the  tips 
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we  expand  ]V}..  'n  ^e  sarne  ^orrn  as  'p*!?,  £  )  given  in  (54).  We 


{here fore  write 


F'<M>  = 


l  Z  «2i.  frr- fi„a> 


Wc  may  expect  this  theory  to  become  inaccurate  for  coses  where  the  influence  of  the 
tips  is  large,  for  example  a  rectangular  plate  of  small  aspect  ratio.  But,  on  the  other 
hand,  small  aspect  ratio  planforms  of  v'ther  shapes  such  as  delta  wings  will  nor  be  in¬ 
fluenced  as  much  by  this  tip  effect,  and  we  therefore  expect  the  approach  to  hold 
even  for  small  aspect  ratios. 

Using  the  solutions  for  pQ(  |  ,  £  )  and  pj(  %  r  ^  )/  eq'->G”or*  (49),  when 
2 

approximated  to  order  ^  ,  becomes 

H  ^  /■»  2 

wtvi5+ty[\  P.(5,p4-+i;%)0‘-S)1p.<5o)‘,5]  +  H[ 


^  f,  11,  j  >4  +  -  (l»|8+V  +  -£ 

-li4)  +^s3"(y"3,al’<-?0-7>P.(i,^4^ 

A  A 

ui.vhii  -5f.«.l)l»c3-?;j,jt-2^U(*.5j. 

A*  A<  "  A% 

P. ».,)!« cj-7?^  - 1 4 \\  «*-s>W )l« J 


«  1  Hl  + 


i t-M'+q-p1-'  Lyv 
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where  G(x~  |  ,  y~  £  )  is  given  in  (46),  In  (59),  the  known  functions  hove  been  co 

lected  on  the  ieft-hond  side  end  oil  the  integrals  there  can  be  integrated  analytically, 
except  the  one  containing  sgn(x  -  ^  ),  which  must  be  integrated  numerically.  The 

integration  of  the  Integra1 
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The  substitution  of  (63)  into  (59)  and  the  use  of  (55)  and  (56)  give 


3 

3TK 


p.iv/nWtt  =  -Lf(r1  + — 2LlL — l&AUx<b 
A  4?t  V  '  (3-  V'  ™ 


(64) 


and 


3Zlt 


[IyL\  +^f  ^5  -  (*"+-£ 

-h  J  -t 


l^4)  ^  ^n(K-5>4(X-5^  >}-r£  )  P»  (5,  £ 

A  A  ■* 

+  -2~rx^\U-5> 


A, 


A 


A, 


+ 

4TC  \  U  1 


(3-  ?> 


:1 

>lJ  typi  l> 


(65) 


where  G(x- ^  ,  y  -  r(  )  is  given  in  (46).  These  two  integral  equations  will  be 

solved  in  the  same  way  as  that  used  to  solve  the  integral  equation  (53).  The  solutions 
may  be  written  as 


M  hi 


f.<m>  -  I  E  <C 


n*  zv  y\ze> 


(66) 


'-•'to 
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end 


M  H 


ft*?  ?)  =  ZL  aZ  in=rJ,u). 


(67) 


The  substitution  c?  (54),  (58),  (66)  and  (67)  into  (63)  gives  us  the  first  four  terms  of  the 
expansion  of  the  pressure  p(  ^  ^  )  up  to  order  ^ Z  c 


5.  Results  for  a  Fiat  Rectangular  Planing  Surface 

We  have  carried  out  numerical  calculations  for  the  case  of  a  flat  rectangular 
plate  for  aspect  ratios,  AR,  varying  from  0.5  to  2.0*.  for  large  aspect  ratios  <he  ex¬ 
pansions  used  in  this  theory  are  not  valid,  while  for  small  aspect  ratios  we  have  rioted 
that  the  theory  becomes  inaccurate  for  planforms  where  the  tip  effect  predominates. 

In  our  presentation  of  the  results  we  use  the  parameter 


-  -fe  '  (68) 

where  2b  is  the  wetted  chord-length.  This  parameter  is  chofifr,  'so  that  the 
resulting  Froucle  number,  1  /j$a  ,  is  based  on  the  characteristic  length  in  the  flow 
direction. 

We  use  the  experimental  data  of  Sottort  (1934)  which  have  sufficient  variation 
in  ^ c  to  make  comparison  with  a  theory  containing  the  gravity  effect  possible. 

In  figure  1  we  show  the  lift  slope  as  a  function  vf  -spect  ratio  for  various  values 
of  P)c  .  The  curve  Ac  =0  represents  half  of  the  lifr  slope  for  the  corresponding 


*  The  aspect  ratio  is  defined  as  the  ratio  of  the  span  to  the  wetted  chord  length  of  the 
planing  surface. 
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wing  wetted  on  both  sides.  Curves  for  ^ c  =0.1  end  0.2  are  chosen  to  show  the  effect 
of  gravity  and  to  correspond  to  SottorPs  values  of  ^ c  .  It  can  be  seen  that  the  theory 
predicts  that  the  effect  of  gravity  increases  the  lift  for  aspect  ratios  less  than  about  1 .5 
and  decreases  the  lift  for  larger  aspect  ratios.  This  tendency  agrees  with  the 
theoretical  results  found  by  Maruo  (1967)  for  the  limiting  cases  of  large  and  smell  as¬ 
pect  ratios  and  with  the  results  of  Sambraus  (1938).  As  the  aspect  ratio  increases  in 
SottorPs  data  the  angle  of  attack  increases  and  decreases.  The  first  grouping 
is  for  around  0.19  with  angle  of  attack  around  4.6°.  The  agreement  between 
theory  and  experiment  is  quite  good.  The  second  grouping  is  for  ^ c  around  C.  16  with 
angle  of  attack  about  5. 1°  and  the  agreement  is  still  good.  As  the  angle  of  attack 
increases  (groups  3,  4  and  5)  the  difference  between  the  theoretical  and  experimental 
results  increases  but  the  agreement  is  fairly  good  overall.  This  indicates  that  the 
non-linear  effect  of  larger  angles  of  attack  may  play  an  important  role.  Generally 
speaking  the  observed  values  are  higher  than  the  calculated  values. 

In  figure  2  we  show  the  position  of  the  center  of  pressure  as  a  function  of 
aspect  ratio  for  various  values  of  |§c  .  The  data  are  much  more  scattered  compared 

with  those  for  the  lift  slope.  However,  agreement  is  reasonable.  Figure  3  illustrates 
the  chordwise  pressure  distribution  at  mid-spon  and  at  90  percent  of  the  span,  for 
aspect  ratios  0.5  and  1 .0.  The  effect  of  gravity  is  to  increase  the  pressure  towards 
the  trailing  edge  and  towards  the  tips.  For  larger  aspect  ratios  the  increase  in  pres¬ 
sure  is  not  as  pronounced. 

In  conclusion,  this  analysis  of  planing  surfaces  of  muderate  aspect  ratio  has 
proved  to  be  reasonably  accurate.  For  smaller  aspect  ratios  we  can  still  obtain  re¬ 
sults  from  this  theory  bur  the  validity  of  the  model  used  is  then  doubtful.  This 
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limitation  is  due  to  the  singular  behavior  cf  the  induced  downwash  near  the  tips,  which 
becomes  increasingly  important  as  the  aspect  ratio  is  reduced.  Further  work  should  be 
directed  towards  the  correct  modelling  of  this  tip  effect. 
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The  steady  motion  of  a  planing  surface  of  moderate  aspect  ratio  at  small  angles  of  attack 
is  considered.  Linearized  theory  is  used  with  a  square-root  type  of  pressure  singularity  re¬ 
presenting  the  flow  near  the  leading  edge.  An  asymptotic  solution  for  the  pressure  distribution 
on  the  planing  surface  at  large  Froude  number  (or  small  B,  the  inverse  of  the  Froude  number)  is 
sought .  The  lowest  order  term  of  the  pressure  distribution,  obtained  by  setting  8  equal  to  zero, 
is  found  to  be  the  same  as  the  pressure  distribution  on  the  lower  side  of  the  corresponding  thin 
wing.  Higher  order  terms  in  Bare  obtained  by  an  iteration  process.  Explicit  solutions  are 
obtained  to  order  B^  for  rectangular  planforms.  Numerical  results  are  calculated  for  rectangular 
flat  plate  planing  surfaces  of  aspect  ratios  from  0.5  to  2.0,  It  is  found  that  for  large  aspect 
ratios  lne  lift  coefficient  is  reduced  by  the  gravity  effect  and  for  small  aspect  ratios  it  is  increased, 
the  dividing  aspect  ratio  being  about  1.5.  The  results  compare  reasonably  well  with  experimental 
data . 
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